Decay dynamics of neutral and charged excitonic complexes in single InAs/GaAs 

quantum dots 
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Systematic time-resoived measurements on neutrai and charged excitonic compiexes (X, XX, X+, 
and XX+) of 26 different single InAs/GaAs quantum dots are reported. The ratios of the decay 
times are discussed in terms of the number of transition channels determined by the excitonic 
fine structure and a specific transition time for each channel. The measured ratio for the neutral 
complexes is 1.7 deviating from the theoretically predicted value of 2. A ratio of 1.5 for the positively 
charged exciton and biexciton decay time is predicted and exactly matched by the measured ratio 
indicating identical specific transition times for the transition channels involved. 

PACS numbers: Valid PACS appear here 



Single quantum dots (QDs) provide the key for quan- 
tum computing [3, 0] and quantum cryptography 0] . For 
real devices, the dynamics of the recombination is of ut- 
most importance, because it limits the maximum mod- 
ulation frequency of the device and the coherence time 
of the emitted light. The electron and hole wave func- 
tions and, therefore, the oscillator strength varies from 
dot to dot because of structural variations Q. Therefore, 
a spread of decay times is expected for different QDs. Ad- 
ditionally, different excitonic complexes in the same QD 
show different decay times for two main reasons: First, 
the number of possible recombination channels varies due 
to differences of the electronic fine structure of the vari- 
ous complexes. Second, the oscillator strength of a spe- 
cific transition depends on the single-particle wave func- 
tions and their overlap being influenced by the number 
of charge carriers constituting the complex. For equal 
specific transition times (STTs), the exciton (X) decay 
time is expected to be twice as long as the biexciton (XX) 
decay time since the XX has two possible decay channels 
and the X has only one. Previous studies of the decay of 
X and XX in InAs QDs were limited to a few QDs and 
reported factors vary from 7.5 [|[ via 2.3 @ to 1.1 0. 
Theoretical predictions range between 4 and 1.5 [1,0. 

In this letter, we present a study of 26 different single 
InAs/GaAs QDs yielding an average factor of 1.7 ± 0.4 
for the decay time ratio of X/XX. Furthermore, we ex- 
pand the study to positively charged excitons (X+) and 
biexcitons (XX+). We predict a ratio of 1.5 and observe 
experimental values yielding 1.5 ± 0.2. This close agree- 
ment shows that the STTs for X+ and XX+ are very 
similar in contrast to the STTs of X and XX. 

The InAs QDs were grown by metal-organic chem- 
ical vapor deposition epitaxy in a GaAs matrix on a 
GaAs(OOl) substrate. During QD growth, the rotation 
of the wafer was interrupted and nominally 1.6 monolay- 
ers of InAs were deposited at 500 °C. The interruption 
leads to a lateral gradient of the QD density. In certain 
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regions of the wafer, the QD density is extremely low 
(«10 7 cm -2 ) giving spectroscopic access to single QDs. 

The QDs were investigated using a JEOL JSM 840 
scanning electron microscope equipped with a cathodo- 
luminescence setup (ioj . The sample was mounted onto 
the tip of a helium flow cryostat providing temperatures 
as low as 6 K. All measurements throughout this letter 
were performed at this temperature. The luminescence 
was dispersed by a 0.3 m monochromator equipped with 
a 1200 lines/mm grating. Single QDs were identified by 
luminescence mapping of the sample. The various QD 
emission lines were a ssig ned to different excitonic com- 
plexes following Ref. [IJJ . 

For time-integrated measurements, a liquid-nitrogen 
cooled Si charge-coupled-device camera was used. The 
time-resolved measurements were performed with a Si 
avalanche photo diode and a beam-blanking unit for 
pulsed excitation. The beam-on time was 5 ns with a 
repetition rate of 25 MHz. This pulse configuration al- 
lowed luminescence decay from a steady-state situation 
thus excluding the influence of the capture process. Per- 
forming transient line shape analysis, the decay times 
could be evaluated with an accuracy of 0.1 ns. The dy- 
namic range of these experiments covered up to three or- 
ders of magnitude enabling the observation of weak slow 
components. 

Typical transients of the decay of four different exci- 
tonic complexes (X, XX, X+, and XX+) in one QD are 
shown in Fig. 1 |l2j. All transients are dominated by a 
fast initial decay typically followed by a second slower 
component. The decay from different QDs shows a sys- 
tematic pattern: The XX+ decays the fastest followed 
by the XX, by the X+, and finally by the X. 

The amplitude of the slower component of the tran- 
sients is very small. For X and X+ its relative amplitude 
compared to the fast decay ranges between and 0.1 and 
for XX and XX+, it ranges between and 0.01. The time 
constants of the slower components vary between 4 and 
9 ns. This component originates from processes feeding 
the initial state of the decay after the external excitation 
has been switched off. X and X+ having longer initial de- 
cay times exhibit a more pronounced second component 
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FIG. 1: Typical transients of the different excitonic complexes 
(X, X+, XX, XX+) of one QD. The transients were fitted by 
a biexponential function. The fast component constitutes the 
radiative decay time. The second component stems from a 
refilling process from neighboring defect states or other QDs. 
The inset shows a typical time-integrated spectrum of a single 
QD. 



than XX and XX+, suggesting a feeding process from an 
external charge carrier reservoir rather than an intrinsic 
feeding. Such slow components were previously contro- 
versially attributed to reemission and lateral transfer of 
charge carriers [l3|, [l4| or to conversion of dark excitons 



to bright excitons via spin- flip processes |15| . For the 
X+, XX, and XX+ decays, the spin flip can be ruled out 
as these complexes possess no such dark states. Instead, 
the source of the feeding process is suggested to be out- 
side of the QDs, e.g., neighboring shallow defect states 
or reemission from other QDs. In the following, we will 
focus on the dominant fast component. 

Figure 2 shows the dominant decay times of X, X+, 
XX, and XX+ as a function of the X+ recombination 
energy. A clear sequence of decay times for the exci- 
tonic complexes is obvious. X exhibits the longest de- 
cay time followed by X+, XX, and finally XX+, analo- 
gous to Fig. 1 and in agreement with previous reports on 
InAs/GaAs QDs [|,[|. 

No clear trend for a dependence of the decay times 
on the X+ energy is identified. Apparently, the struc- 
tural parameters, which lead to differences in the tran- 
sition energies, do not control the decay times in the 
same way. Interestingly, the scatter of the X decay times 
is much larger than that of the other excitonic com- 
plexes. The mean decay times are t(X) — 1.22 ±0.25 ns, 
t(X+) = 0.97 ± 0.15 ns, t(XX) = 0.76 ± 0.12 ns, and 
t(XX+) = 0.66 ± 0.14 ns. The standard deviation of 
t(X) is almost twice as large as for the other complexes. 
The X wave function seems to be more sensitive to the 
structural properties of the QD than the wave functions 
of complexes containing more holes. 

The decay times of the neutral and charged complexes 
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FIG. 2: The dominant decay times of the different excitonic 
complexes are shown as a function of the X+ recombination 
energy. No clear dependence on energy can be seen. 



will now be compared with respect to the number of al- 
lowed transitions of each excitonic complex. 

The radiative decay time t(EC) of an excitonic com- 
plex EC is composed of the specific transition time Tj / of 
all recombination channels involved flql: 
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with the mean population m (X) n * = 1) of the initial 
state. It is assumed that the STTs for different tran- 
sition channels — > \f) and — ► |/'}) of one com- 
plex are identical (jij — iv,// = tec) since the oscilla- 
tor strength is independent of the actual spin configura- 
tion as long as the spin selection rules are obeyed. Un- 
der this precondition, the expected excitonic decay times 
t(X), t{X+), t(XX), and t(XX+) are determined by 
the number of allowed recombination channels and their 
STTs, respectively 

Figure 3 shows the energy schemes of the four excitonic 
complexes to illustrate the different transition channels. 
The initial and final states are labeled with their respec- 
tive total spin configurations. 



1-1>-|1> 
l-D+ID 

i-2>-|2) 
|-2>+|2> 



L +1,0 



XX 



x+ 



l 0,+1 



|0> 



"0 : -1 



I-1>-I1> 



|±1/2> 



XX+ 



■|0> 



|-1>+|1> 

J;2>-|2> 
T-2>+|2> 



T +1/2, 
+3/2 

T -1/2, 
-3/2 



T +3/2, 
+1/2 

T -3/2, 
-1/2 



|±3/2> 



— |±3/2) 

T +3/2, 
+5/2 

T -3/2, 
-5/2 

— 1±1/2> 
-|±5/2> 
_|±7/2> 



FIG. 3: Energy scheme for the recombination channels of 
the different excitonic complexes. The spin configurations of 
the different initial and final states are shown. Thick lines 
correspond to twofold degenerate states. 
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FIG. 4: Ratio of the decay times of the neutral and charged 
complexes as a function of X+ recombination energy. The 
dashed lines indicate the mean values. The means of the ratios 
are 1.7±0.4 and 1.5±0.2 respectively. The greater scatter 
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ratio stems from the greater scatter of the X 



The initial configuration of the optically active X con- 
sists of the two bright states with the spin configuration 

| - 1) ± |1) 0, m The dark states I - 2 ) ± I 2 ) are 
optically forbidden and do not contribute to the recom- 
bination process. The exciton fine structure splitting of 
the QDs probed here was determined to be less than 
20 fj,eV. Hence, the mean initial population n±± of the 
two bright states is \ . Each state possesses one recombi- 
nation channel to the ground state |0). The decay time 
then amounts to t(X) = Tx- 

In contrast to X, XX consists of a single initial state 
|0) with no = 1 and possesses two recombination chan- 
nels to the two final states: the bright states of the 
X (| — 1) ± |1)). Hence, the decay time amounts to 
t(XX) = \t xx . 

To compare the decay times of X and XX, identical 
STTs for both complexes are assumed. Narvaez et al. 
fil| calculated almost identical STTs (tx = txx) f° r 
the X and XX in lens shaped InAs/GaAs QDs of vari- 
ous sizes. Under this assumption, a decay time ratio of 
t Syy\ — 2 follows. 

The initial and final states of X+ consist of two de- 
generate states | ± i) and | ± |), respectively. For each 
initial state, only one final state is allowed due to the 
spin selection rules. Since the initial states are equally 
populated (n±i/2 = o)> the decay time of X+ following 
Eq. (|T|) yields: 
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The initial configuration of XX+ consists of the two 
degenerate states | ± |). The final state of the XX + 
is a triplet, each component consisting of three doublets 
[ljj[2(|. The transition to the | ±7/2) state is forbidden, 



while the transition to the | ± 5/2) state is allowed. The 
transition probability into the | ± 1/2} state is reduced 
to \ since the final state includes a forbidden transition 
[2l| . Thus, the decay time of the XX+ accounts to 
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Just like for the neutral complexes, equal STTs (tx+ = 
T xx+) for the charged complexes are assumed. Since the 
probed QDs are of small size and, therefore, in the strong 
confinement regime, adding a neutral pair of charge car- 
riers should alter the wave functions only slightly. The 
expected ratio of decay times J(xx+) then amounts to 
1.5. 

In Fig. 4, the measured ratios of the X and XX de- 
cay times (black squares) and of the X+ and XX+ decay 
times (red circles) are plotted as a function of X+ recom- 
bination energy. Again, no apparent trend can be seen. 
The mean value of J(x~x) is 1-^ with a standard deviation 



of 0.4. For 



r(X±) 



the mean value results to 1.5 with a 



r(XX + ) 

standard deviation of 0.2. 

The mean J^ X x) vam e is clearly below the predicted 
value of 2 with a large scatter of the individual values. 
Qualitatively, the argument based on the number of in- 
volved transition channels holds. The exact value, how- 
ever, crucially depends on the specific QD since the STTs 
differ a lot from case to case and are not always compara- 
ble as in [16(] . Wimmer et al. calculated a size dependent 
ratio between 1.5 and 2 [8]. The change in the ratio is 
governed mainly by t(X), thus, supporting our measure- 
ments. 

The measured mean value of J^xx+) matches exactly 
the expected value of 1.5. The scatter is considerably 
smaller than for T ,i?2-, with a standard deviation of 0.2 
indicating a stronger similarity of the STTs of X+ and 
XX+ from dot to dot than for the neutral complexes. 
Here, the number of allowed transitions is the key pa- 
rameter determining the ratio of decay times of both 
complexes. 

The comparison of the measured t(X) and t(X+) 
shows that the STTs tx and tx+ differ strongly. Adding 
a single charge carrier to an excitonic complex influences 
the wave functions and, therefore, their overlap is much 
stronger than adding a neutral pair of charge carriers. 

In summary, we presented a systematic study of the de- 
cay dynamics of excitonic complexes in 26 different single 
InAs/GaAs QDs. The influence of the specific transi- 
tion times and electronic fine structure on the dynamics 
was analyzed. Surprisingly, the scatter of the neutral 
exciton's decay time from dot to dot is larger than the 
scatter of the other complexes' decay times indicating a 
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greater sensitivity of the exciton wave function to details 
of the structural properties of the QD. No dependence on 
the transition energy was observed in the 1.26 - 1.35 eV 
range. The mean value of J^xx) was ^ thus, 
below the value of 2 as given by the different number of 
decay channels alone. This implies shorter STTs for X 
than for XX. For J^xx+) a vam e of 1.5 was predicted. 
The mean measured ratio is 1.5 ± 0.2, thus, matching the 
expected theoretical value. This indicates similar STTs 



for the transition channels of X+ and XX+. The addi- 
tional positive charge carrier seems to stabilize the wave 
function overlap when an additional exciton is added to 
the complex. 

This work was supported by the Deutsche Forschungs- 
gemeinschaft in the framework of SfB 787 and the 
SANDiE Network of Excellence of the European Com- 
mission, Contract No. NMP4-CT-2004-500101. We are 
indebted to K. Potschke for the growth of the sample. 



[1] G. Chen, T. H. Stievater, E. T. Batteh, X. Li, D. G. 

Steel, D. Gammon, D. S. Katzer, D. Park, and L. J. 

Sham, Physical Review Letters 88, 117901 (2002). 
[2] P. Bianucci, A. Muller, C. K. Shih, Q. Q. Wang, Q. K. 

Xue, and C. Piermarocchi, Physical Review B 69, 161303 

(2004). 

[3] O. Benson, C. Santori, M. Pelton, and Y. Yamamoto, 
Physical Review Letters 84, 2513 (2000). 

[4] A. Schliwa, M. Winkelnkemper, and D. Bimberg, Physi- 
cal Review B 76, 205324 (2007). 

[5] M. H. Baier, A. Malko, E. Pelucchi, D. Y. Oberli, and 
E. Kapon, Physical Review B 73, 205321 (2006). 

[6] R. M. Thompson, R. M. Stevenson, A. J. Shields, I. Far- 
rer, C. J. Lobo, D. A. Ritchie, M. L. Leadbeater, and 
M. Pepper, Physical Review B 64, 201302 (2001). 

[7] C. Zinoni, B. Alloing, C. Monat, V. Zwiller, L. H. Li, 
A. Fiore, L. Lunghi, A. Gerardino, H. de Riedmatten, 
H. Zbinden, and N. Gisin, Applied Physics Letters 88, 
131102 (2006). 

[8] M. Wimmer, S. V. Nair, and J. Shumway, Physical Re- 
view B 73, 165305 (2006). 
[9] G. A. Narvaez, G. Bester, A. Franceschetti, and 
A. Zunger, Physical Review B 74, 205422 (2006). 
[10] J. Christen, M. Grundmann, and D. Bimberg, Journal of 

Vacuum Science and Technology 9, 2358 (1991). 
[11] S. Rodt, A. Schliwa, K. Potschke, F. Guffarth, and 

D. Bimberg, Physical Review B 71, 155325 (2005). 
[12] The XX+ transient was always obtained from the more 



intense of the two emission lines. The other line yields 
the same decay time since the occupation of the same 
complex is probed. 

[13] S. Rodt, V. Tiirck, R. Heitz, F. Guffarth, R. Engelhardt, 
U. W. Pohl, M. Strassburg, M. Dworzak, A. Hoffmann, 
and D. Bimberg, Physical Review B 67, 235327 (2003). 

[14] A. S. Shkolnik, L. Y. Karachinsky, N. Y. Gordeev, G. G. 
Zegrya, V. P. Evtikhiev, S. Pellegrini, and G. S. Buller, 
Applied Physics Letters 86, 211112 (2005). 

[15] J. M. Smith, P. A. Dalgarno, R. J. Warburton, A. O. 
Govorov, K. Karrai, B. D. Gerardot, and P. M. Petroff, 
Physical Review Letters 94, 197402 (2005). 

[16] G. A. Narvaez, G. Bester, and A. Zunger, Physical Re- 
view B 72, 245318 (2005). 

[17] M. Bayer, G. Ortner, O. Stern, A. Kuther, A. A. Gor- 
bunov, A. Forchel, P. Hawrylak, S. Fafard, K. Hinzer, 
T. L. Reinecke, S. N. Walck, J. P. Reithmaier, 
F. Klopf,and F. Schafer, Physical Review B 65, 195315 
(2002). 

[18] The normalization factor \[2 has been disregarded 

throughout this letter for simplicity. 
[19] R. Seguin, S. Rodt, A. Schliwa, K. Potschke, U. W. Pohl, 

and D. Bimberg, phys. stat. sol. (b) 243, 3937 (2006). 
[20] K. V. Kavokin, physica status solidi (a) 195, 592 (2003). 
[21] I. A. Akimov, K. V. Kavokin, A. Hundt, and F. Hen- 

neberger, Physical Review B 71, 075326 (2005). 



